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The origin of lithium (Li) and its production process have long been an unsettled question 
in cosmology and astrophysics. Candidates environments of Li production events or sites 
suggested hy previous studies include hig hang nucleosynthesis, interactions of energetic cos¬ 
mic rays with interstellar matter, evolved low mass stars, novae, and supernova explosions. 
Chemical evolution models and observed stellar Li abundances suggest that at least half of 
the present Li abundance may have been produced in red giants, asymptotic giant branch 
(AGB) stars, and novae However, no direct evidence for the supply of Li from stellar ob- 
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jects to the Galactic medium has yet been found. Here we report on the detection of highly 
hlue-shifted resonance lines of the singly ionized radioactive isotope of heryllium, ^Be, in the 
near ultraviolet (UV) spectra of the classical nova V339 Del (Nova Delphini 2013). Spectra 
were obtained 38 to 48 days after the explosion. ^Be decays to form ^Li within a short time 
(half-life 53.22 days"). The spectroscopic detection of this fragile isotope implies that it has 
been created during the nova explosion via the reaction ^He(a,7)^Be (ref. 5), and supports 
the theoretical prediction that a significant amount of ^Li could be produced in classical nova 
explosions. This finding opens a new way to explore ^Li production in classical novae and 
provides a clue to the mystery of the Galactic evolution of lithium. 

V339 Del (= Nova Delphini 2013) is a elassieal nova that was diseovered as a bright 6.8 
magnitude (unfiltered) souree on 2013 August 14.584 Universal Time (UT)^. After 40 hours from 
the diseovery, a maximum was reaehed on Aug 16.25 (MJD = 56,520.25) at U = 4.3 (ref. 7). Then, 
it began a normal deeline. 

High-resolution speetra (R = 90, 000-60, 000) of V339 Del were obtained at four epoehs 
after its outburst (-1-38, -i-47, -i-48, and -i-52 d). These speetra eontain a series of broad emission 
lines originating from neutral hydrogen (H l, Balmer series) and other permitted transitions of 
neutral or singly ionized speeies (e.g.. Fell, He I, Call). These emission lines are usually seen 
in post-outburst speetra of elassieal novae. Most of these broad emission lines are aeeompanied 
by sharp and blue-shifted multiple absorption lines at their blue edges. The typieal radial veloeity 
(Uad) of these highly blue-shifted absorption lines is ~ —1, 000 kms“^. Figure 1-a and -b display 
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the spectrum obtained at +47 d in the vicinities of H and Ca ll K lines, respectively. The H l line 
is accompanied by a broad emission with a FWHM of ~ 1,300 km s“^ centered at Vrad ~ 0 km s“^. 
The Ca ll K (and H) has a weak but broad emission and strong absorption components caused by 
the interstellar (IS) absorptions. In addition to these profiles around the rest positions of both lines, 
two sharp absorption components are found at Vrad = —1, 268±2 and —1,103±1 kms“^ Among 
these, the latter is apparently stronger than the former. Such absorption line systems have been 
found in post-outburst spectra of several classical novaeThe absorption line systems in V339 
Del contain numerous transitions originating from singly ionized iron-group species (Fe ll, Ti It, 
Cr II, Mn II, and Ni II). The depths of all blue-shifted absorption lines in V339 Del are only ^ 25% 
of the continuum, while the bottoms of some strong lines (e.g., Balmer lines; see Extended Data 
Fig. 2) have flat features suggesting that the absorption is saturated. These observational results can 
be interpreted as the effect of absorbing materials partially covering the background light source. 
There are no Na l D doublet lines, which are often found to be the strongest absorption features in 
novae within a few weeks after their outbursts^’'*’. We interpret this as indicating that the ionization 
state of the ejected gas has evolved into a higher stage of excitation before our observing epochs 
(5-7 weeks). The observed spectral energy distribution of this nova indicates that the shape of 
the continuous radiation had entered a very hot stage (effective temperature >100,000 K) within 5 
weeks after the explosion Other observed characteristics of this nova (e.g., light curves, optical 
and UV emission lines) show that it is a typical Fe ll nova with a CO white dwarf (WD) 

Among these absorption line systems, we have noticed two remarkable pairs of absorption 
features near 312 nm. These correspond to the absorption components originating from transitions 
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at ~313 nm. These pairs are marked as A, B and C, D, respeetively, in Figure 1-e. Adopting the 
wavelengths of the resonanee doublet lines of singly ionized ^Be at 313.0583 and 313.1228 nm 
we find that features A and B eoineide with the —1,103 kms“^ eomponents of the 313.0583 and 
313.1228 nm lines, respeetively. Similarly, features C and D eoineide with the —1,268 km s~^ 
eomponents of these two lines. Separations in wavelength between features A and B, and between 
features C and D are eonsistent with the separation between the doublet lines within the measure¬ 
ment uneertainties. Figure 1-d illustrates these eoineidenees on the veloeity seale. Thanks to the 
high resolution of the speetrum (~ 0.0052 nm), we ean elearly distinguish them from the doublet 
of ®Be II at 313.0422 and 313.1067 nm'"^. After ruling out the possibilities of alternative identifi- 
eations, we eonelude that these absorption features at 312 nm are eaused by ^Be, and not by ®Be 
whieh is the only stable isotope of Be. Original ®Be eontained in the progenitor star would have 
been depleted during its evolution beeause this isotope is destroyed at temperatures T > 3 x 10® 
K. On the other hand, produetion of the unstable isotope ^Be by the reaetion ^He(Q;, 7 )^Be in nova 
explosions has been theoretieally predieted 

The transition probability of the ^Bell line at 313.0583 nm {loggf = —0.178) is twiee as 
large as that of the ^Be ll at 313.1228 nm (log gf = —0.479) Due to saturation effeets, the ratio 
of their equivalent widths is expeeted to be in the range between 2 (no saturation) to 1 (eomplete 
saturation). The measured ratios are 1.1 ± 0.3 and 1.6 ± 0.4 for the eomponents at Wrad = —1, 268 
and —1,103 km s“\ respeetively. These are within the range expeeted for the doublet, although the 
values eontain some errors (^ ±25%) due mainly to the uneertainty in the eontinuum plaeement. 
The weaker eomponent at Urad = — 1, 268 km s“^ has a ratio eloser to eomplete saturation. This ean 
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be interpreted as resulting from the fact that the absorbing gas cloud moving with Urad = — 1, 268 
km s“^ has a smaller covering factor, and at the same time, higher column density of ^Be ion than 
the gas cloud moving with Urad = —1,103 km s~^ 

Figure 2 displays the velocity plots of normalized spectra for different species of the ab¬ 
sorption line systems at four observing epochs from -i-38 to -i-52 d. On -i-38 d, the ^Be ii doublet 
has an absorption component at = —1, 386 ± 3 kms“^and shows a complicated profile near 
^'rad ~ —1,000 kms“^. From +A1 to -i-48 d, the absorption components at Wrad = —1,268 and 
—1,103 kms“^ in -i-47 d shift by —26 ± 3 and —17 ± 4 kms“^ blueward, respectively. These 
changes can be interpreted as due to the fact that we are observing accelerating blobs of nova 
ejecta. All of the blue-shifted absorption line systems had disappeared in the spectrum of -i- 52 
d except for the meta-stable Hel lines at 318.8 and 388.9 nm. This fact indicates that the gas in 
the absorption line systems have evolved into a higher ionization stage as discussed in the case of 
the nova VI280 Sco^^. These observations show - (1) Several blue-shifted absorption lines with 
different velocities are found from different species at each epoch; (2) Radial velocities of differ¬ 
ent transitions belonging to a velocity component determined by Gaussian fittings agree within 
lAwradl <1-3 kms“^; (3) Each component shifts blueward with time, indicating that the ejecta is 
being accelerated; (4) The strengths of the blue-shifted absorption lines weaken quickly during the 
observing period. The velocities and the strengths of the ^Be ll doublet behave perfectly synchro¬ 
nized with those of other species. This means that the gas producing the absorption line systems of 
V339 Del contains a considerable amount of ^Be ion which can produce detectable absorption lines 
and strongly suggests that the gas must have experienced an explosive thermonuclear runaway on 
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the surface of the WD. 


Our spectroscopic detection of ^Be in a classical nova immediately connects to the produc¬ 
tion of ^Li. The production of ^Be via the nuclear reaction ^He(Q;, 7 )^Be in novae has been studied 
theoretically. However, no observational confirmation has been made up to the present time. This 
is because the presence of ^Be is very transient and is only observable in the near UV range where 
the atmospheric absorption severely obstructs observations from ground-based telescopes. In the 
case of V339 Del, the ^Be doublet can be identified only within a very short period (from ~6 to 
~7 weeks after the outburst). In earlier epochs, it might be difficult to identify due to saturation ef¬ 
fects. For nearby bright novae, there have been several attempts to detect the 478 keV gamma-ray 
line produced by the decay of ^Be. However, no definite detection has been reported because of 
the insufficient sensitivity The ^Be absorption lines in the near UV spectrum of V339 Del are 
found in highly blue-shifted (~ 1, 000 kms“^) flows, which have been blown off by the outburst. 
This means that it will soon decay to ^Li in cooler interstellar or circumstellar matter on a time 
scale given by the half-life of ^Be (53.22 days). The absence of the ^Lil line at 670.8 nm in our 
spectra can be interpreted as due to the fact that all of Li in the absorbing material of V339 Del has 
been ionized during the observing period as mentioned above. This is in accordance with the fact 
that no Na l D lines are found in the absorption line system. 

The ^Be ii doublet corresponds to the doublet of the Ca ii resonance lines (H and K lines) on 
the atomic energy level diagrams. The CallK line at 393.366 nm has \og gf = -1-0.135 (ref. 14). 
Supposing that most of the ^Be and Ca in the absorption line system are singly ionized and the 
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resonance lines of both ions are unsaturated, the ratio of their equivalent widths directly reflects 
the number density ratio between ^Be and Ca ions. In the spectrum obtained at +47 d, the ratios 
of ^Bell (313.1228 nm)/CallK, which are less affected by saturation and/or contamination, are 
~ 1.3 ± 0.3 and ~ 0.7 ± 0.2 for the —1, 268 and the —1,103 kms“^ components, respectively. 
This means that the column number density of ^Be is 5.3-2.9 times higher than that of Ca. Using 
this ^Be/Ca number ratio, the mass fraction of ^Be relative to the sum of all the constituent mass 
components, X(^Be), is presented as (4.4 ± 2.2) x 7/40 x X(Ca). If we adopt the solar X(Ca) 
(= 10“"^'^®), for instance, the X(^Be) in the absorbing gas system should amount to ~ 

The error estimation includes the uncertainty in the local continuum placement (^ ±25%) and the 
difference of derived equivalent widths derived from individual velocity components (~ ±30%). 
Several additional factors, such as the abundance of Ca in the absorbing gas, and the difference in 
the ionization state between ^Be and Ca, are difficult to estimate because the nova ejecta model is 
not yet established. Taking such uncertainties into account, the error involved in the above estimate 
of the ^Be abundance could be even larger by a factor of several. However, in spite of the large 
uncertainty, the abundance of ^Be is larger than, or at least as large as, theoretical predictions for 
CO novae [e.g., X(^Be) ^ 10“^-^ (ref. 21)]. This indicates that classical novae could play an 
important role as contributors of ^Li in the Galaxy. 

The observed ^Li evolutionary curve^ has a plateau for young Galactic ages (^ 2.5 Gyr) 
followed by a steep rise. To explain this requires a relatively low-mass stellar component that 
evolves over a long lifetimes. Candidates for this, such as low-mass red giants or novae, have 
been proposed to be major sources of ^Li production (> 50% of the solar system Li measured 
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in meteorites) in the Galaxy The production of ^Li in low-mass stars has been theoretically 
studied and Li-enhanced red giants and AGB stars are indeed identified The contribution 
to Li enrichment in the Galaxy by these objects has, however, not been confirmed. The reason for 
this is that the Li-rich phase in these stars might be of quite limited duration and the contribution 
is dependent upon the mass-loss rate of such objects. Nova eruptions involve a long delay time 
before working as stellar ^Li factories. This is because ^He rich low mass secondaries are necessary 
to produce ^Be efficiently via the ^He(Q;, yj^Be reaction^*. It is important to know whether this 
phenomenon is common among classical novae to quantify their contribution to the rapid increase 
of ^Li in the Galaxy. Since V339 Del appears to be one of the ordinary Fell type novae which 
occupy ~60 % of all classical novae the ^Be production found in this object might be occurring 
in many classical novae. Our successful detection of ^Be in V339 Del indicates that measurements 
of the ^Be lines in the near UV range for post-outburst novae within the lifetime of this isotope is 
a powerful way to estimate the contribution of novae to the chemical evolution of lithium in the 
Galaxy. 
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Figure 1 : Blue-shifted absorption line systems in the spectrum of V339 Del obtained at +47 

day. Panels a-c display the speetrum in the vieinity of H r/ (a), Ca ii K (b), and the ^Be ii doublet 

(c), on the veloeity (upper horizontal) scale. Two blue-shifted absorption components and the 

zero velocity position for each line are indicated by ticks, c, The velocity scale is adjusted to one 

of the ^Bell doublet (313.1228 nm, red). The positions of blue-shifted components of Cm at 

313.205 nm are displayed at the bottom. Panel d shows the enlarged radial velocity profiles. The 
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vertical dashed lines show two common absorption components. The expected positions of the 
®Be II doublet are indicated by green lines. 
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Figure 2: Time variations of the blue-shifted absorption line systems from +38 to +52 d. 

Absorption line systems originating from different species at four observing epochs are plotted 
on the velocity scale. All lines are normalized to the local continuum. Blue-shifted absorption 
components observed at each epoch are indicated by vertical dashed lines. The identified line or 
expected line contaminations are labeled above each lines, a, On +38 d, the profile of the ^Beii 
doublet around v^sa ~ —1, 000 kms“^ is complicated, and possibly interpreted as being saturated, 
d, On +52 d, no blue-shifted absorption can be found except for the metastable He l lines. 
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Methods 


Discovery V339 Del (= Nova Delphini 2013) is a classical nova that was discovered as a bright 
6.8 magnitude (unfiltered) source by Koichi Itagaki on 2013 August 14.584 UT and announced 
in American Association of Variable Star Observers (AAVSO) Alert Notice^. Its progenitor is 
estimated to be a blue star USNO B-1 1107-0509795 with B ~ 17.20, Rc ~ 17.45 on the first 
Palomar Sky Survey Plates (exposed on 1951 July 7), and with B ~ 17.39, Rq ~ 17.74 on the 
second Palomar Sky Survey plates (exposed on 1990 July 18 and September 15, respectively) 
No significant changes were found in its photometric behavior for at least a few years prior to the 
outburst^^. On an unfiltered pre-discovery image obtained on 2013 August 13.998 UT, the object 
was still at 17.1 mag^^. This means that the object was still in quiescence until at least 14 hours 
before its discovery, and that it showed a very fast rise to the maximum. After 40 hours from the 
discovery, maximum was reached on Aug 16.25 (MJD = 56,520.25) at V = 4.3 (ref. 7). Then, it 
began a normal decline. The nova had been detected as a transient high energy gamma ray (>100 
MeV) source within ~10 days after the outburstAngular sizes of the expanding shell around the 
nova had been monitored until ~-i-40 d using near-infrared interferometric observations^^. Then, 
combining the expansion velocity obtained in the optical region, the distance to the nova had been 
derived as 4.54 ± 0.59 kiloparsecs from the Sun. 

Observations and data reduction Post-outburst spectra of V339 Del were obtained using the 
High Dispersion Spectrograph (HDS)^^ of the 8.2 m Subaru Telescope at four epochs from 2013 
September 23 to October 7 (-1-38, -i-47, -i-48, and -i-52 d after the maximum). According to the 
AAVSO light curves (see in Extended Data Fig. 1), our first observation was just before the start 
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of the rapid decline in optical magnitudes by dust formation The following three were obtained 
during the continuous decline. We obtained spectra under 3 configurations of the spectrograph, 
which cover the wavelength regions from 303 to 463 nm, from 411 to 686 nm, and from 667 to 
936 nm. Spectral resolving power was set to i? ~ 90, 000 or 60,000 with 0".4 (0.2 mm) or 0".6 
(0.3 mm) slit widths, respectively. The exact times and wavelength ranges of obtained spectra are 
summarized in Extended Data Table 1. Data reduction was carried out using the IRAF software 
in a standard manner. The non-linearity in the detectors are corrected by the method described in 
ref. 38. The wavelength calibration has been performed using a Th-Ar comparison spectrum and 
the typical residual in wavelength calibration is ^ 10“"^ nm (~0.1 kms“^) for each spectrograph 
configuration. The typical systemic variance of the spectrograph is ^ 10“^ nm per an hour. We 
also examined the accuracy of radial velocity determination in our measurement using the iden¬ 
tified iron-group transitions in 315-351 nm. For the spectrum obtained at -i-38d, the velocity of 
the strongest component in the absorption line system was —996.1 ± 0.7 kms“^ determined by 
Gaussian fittings. In total, we concluded that the residual in our velocity scale determination was 
~ ±1 kms“^. Spectrophotometric calibration was performed using the spectrum of BD-i-28° 4211 
(ref. 39) obtained nearly at the same altitude of the nova on the same nights. All spectra were 
converted to the heliocentric scale. Correction for interstellar extinction has not been applied. The 
average signal-to-noise ratio in the spectra obtained at four epochs is ~140 at ~312 nm, where we 
found the ^Be lines. 

Highly blue-shifted absorption line system The spectra of V339 Del exhibit a series of broad 
Fe II emission lines, which indicate that the object is a typical Fe ll type nova^^. Since no strong 
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emission originating from Ne is found in the spectrum even in +52 d, the WD in the system is 
supposed to be a CO-WD 

Extended Data Fig. 2-a displays the radial velocity profiles of three Fe ll belonging to multi- 
plet number^® (42) in the spectrum of +38 d. The absorption line system on +38 d clearly consists 
of five components at Wrad = —954, —968, —985, —996, and —1, 043 kms“^ as indicated by the 
dashed lines in Extended Data Fig. 2-b. Similarly, blue-shifted absorption components are found in 
Baimer lines (Extended Data Fig. 2-c) and also in other permitted lines (Ca II H and K, He I at 587.6 
nm). In the near UV range, numerous absorption lines in the complex continuum are identified as 
the transitions of singly ionized iron-group species. Most of them belong to the absorption line 
systems found in the visual region (Extended Data Fig. 3). We applied a Doppler correction using 
the the radial velocity of the strongest blue-shifted absorption line in the system to identify sources 
of transitions (see in Extended Data Fig. 3-b). We use the velocities for Doppler corrections as 
V+S8 = —996, V+4Y = —1,103, and = —1,120 kms“^ for +38, +47, and +48 d, respectively 
(in Extended Data Fig. 3-b and Fig. 4-a, b, and c). All of the identified transitions originate from 
levels of low excitation potentials (^4 eV). The residual intensity at the bottom of these lines ex¬ 
ceeds 75 % of the continuum, while the bottoms of some strong Fe ll and Balmer lines show flat 
structures. These observations suggest that the saturated absorption lines are created by clouds of 
absorbing gas, which cover the continuum emitting region by only about 25 %. 

Very similar short-lived blue-shifted metallic absorption systems (Transient Heavy Element 
Absorption; THEA) have been reported in post-outburst spectra of several classical novaeEs¬ 
pecially, great majority of novae show strong blue-shifted (400-1,000 kms“^) multiple absorption 
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components in the NalD doublet in days following their outbursts. In the very slow nova VI280 
Sco, multiple high-velocity (700-900 kms“^) absorption components have been found in the Nal 
D doublet even 800 days after its maximum^. Although no absorption components of the Nal D 
are found at all epochs of our observations of V339 Del, other characteristics of the absorption 
line systems found in V339 Del are quite similar to those of the THEAs in other novae. They are 
(a) highly blue-shifted (~ 1, 000 kms“^), (b) divided into several velocity components, (c) time 
variable in their shapes and velocities, and (d) short-lived (2-8 weeks). 

Contamination to the ^Be ii doublet We carefully inspected possible contaminations of absorp¬ 
tion lines originating from other species to the ^Bell lines consulting the atomic line database"^^ 
Extended Data Eig. 4 displays the spectra in the vicinity of the ^Be ll doublet obtained at four 
epochs of our observations. 

On the spectrum obtained at -i-47 d, there are no candidates to contaminate the —1,103 or 
the —1,268 kms“^ components of ^Bell at 313.1228 nm, which we use in our ^Be abundance 
estimation (see in Extended Data Eig. 4-b). 

At this epoch, the other line of the ^Be doublet at 313.0583 nm may be contaminated by 
some lines originating from iron-group species. We estimate that the —1,268 kms“^ component 
of Cm (5) at 313.205 nm {loggf = -1-0.079) may contaminate to the —1,103 kms“^ component 
of this ^Be ll line. The influence of this contamination can be evaluated adopting the line strength 
ratio between the pair of velocity components of Cm (5) at 312.497 nm (log gf = -1-0.018) to that 
of Cm (5) at 313.205 nm. It is quite small compared with the strength of the ^Be ll line (^ 5 %). 
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Concerning the —1, 268 kms“^ component of this ^Be ll line, we conclude that the weak lines of 
Fell(96) at 312.901 nm (loggf = —2.70) and Cm(5) at 312.869 nm (loggf = —0.32) do not 
contaminate severely. This is because similar weak lines of Fell(82) at 313.536 nm (loggf = 
—1.13) and Cm (5) at 313.668 nm (log gf = —0.25) had completely disappeared until + 47 d. We 
can neglect the contamination from the V ll(l) line at 313.0257 nm (log gf = —0.29), because the 
other Vll(l) line at 312.621nm (loggf = —0.27) is not detected on our spectra. 

^Be abundance estimation We empirically estimate the abundance of ^Be in the absorbing gas 
by comparing the equivalent widths of the ^Be ll line with those of the Ca ll K line that are the 
similar transitions on the atomic energy level diagrams. We assume that the covering factor of the 
absorbing gas cloud to the background illuminating source has no wavelength dependence. This 
method could be a simple and robust approach to estimate the abundance ratio independent of 
ejecta models for nova explosions. The estimate, however, includes some uncertainties. One is 
the difference in the ionization potentials between Be (the first and second ionization potentials; 

11 = 9.32, I 2 = 18.21 eV) and Ca (Ji = 6.11, I 2 = 11.87 eV)'"^ that could result in a difference 
of ionization states between Be and Ca. However, all of iron-peak elements (Ti to Fe) found 
in the absorption line systems, which have intermediate ionization potentials (Ji = 6.75-7.90, 

1 2 = 13.58-18.12 eV) between those of Be and Ca, are observed only in singly ionized states, 
suggesting that dominant fractions of Be and Ca are in the singly ionized states, too. In the obtained 
spectra, we could not find any resonance lines of Sril or Ball, which correspond to those of 
Bell and Can. The Sr/Ca and the Ba/Ca number ratios would be quite small, as seen in the 
solar abundance (<^ 10“^). Another uncertainty is the X(Ca) in the absorption line system. Our 
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assumption that the absorbing gas has the solar X (Ca) would not be far from the reality, because 
the theoretical analysis predicts no overabundance of elements with the mass number >30 in ejecta 
of CO novae^^ We remark that our ^Be abundance estimation is carried out using the data obtained 
at +47 d, which is close to the half-life of ^Be (53.22 days). Therefore, the abundance of the freshly 
produced ^Li in this nova explosion could be ~2 times higher than the X (^Be) on +47 d. 
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Extended Data Table 1: Journal of HDS observations of V339 Del 


Extended Data Eig. 1: Optical light curves of V339 Del. V (green) and R (red) magnitudes are 
taken from the AAVSO database. The epochs of our HDS observations are indicated by arrows. 


Extended Data Eig. 2: Optical spectrum of V339 Del obtained at +38 d. Panel a shows the 
radial velocity plots of three Ee ll emission lines belonging to the same multiplet number^® (42). 
In addition to the similarity of their broad emission profiles, all lines have common blue-shifted 
absorption line features around —1, 000 kms“^. Panel b shows the enlarged view of the 

absorption line features in panel a. Dips of individual absorption line are indicated with dashed 
lines, c, The absorption line systems in Hi Balmer lines drawn on the same velocity scale as in 
panel b. 


Extended Data Eig. 3: Near UV spectrum of V339 Del obtained at +38 d. Panel a shows the 
overall view of the spectrum from 308 to 350 nm. Identified Ee ll emission lines are indicated with 
red ticks at the bottom. The identified absorption line systems originating from iron-group ions - 
Eell (red), Tin (blue). Cm (green), Mnll, Nill, and Vll (black) - are indicated by ticks at the 
top. Panel b shows a sample of the absorption line identification. The results of our identification 
are displayed along the spectrum, c. Same as Extended Data Eig. 2-b, but for two lines (Tin and 
Cr II) highlighted in panel b are plotted on the velocity scale. 
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Extended Data Fig. 4: Spectra in the vicinity of the Be ii doublet from +38 to +52 d. a - c, The 

horizontal scale is displayed with the heliocentric (bottom) and the Doppler corrected wavelengths 
(top). The Doppler corrections are applied using Wdays = n+ 47 , and n +48 for panels a, h, and c, 

respectively. The local continuums fitted with high order (10-20) spline functions are over-plotted 
with green lines. The positions of the strongest (Vrad = todays) and the second strongest components 
of the absorption line system are indicated by colored long and short lines connected by horizontal 
bars, d, Since no apparent absorption lines are found in +52 d, the spectrum is applied a Doppler 
correction using n+ 48 . 
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Extended Data Fig. 1: Optical light curves of V339 Del. 
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Extended Data Fig. 2: Optical spectrum of V339 Del obtained at +38 d. 
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Extended Data Fig. 3: Near UV spectrum of V339 Del obtained at +38 d. 
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Extended Data Fig. 4: Spectra in the vicinity of the Be ii doublet from +38 to +52 d. 
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Extended Data Table 1: Journal of HDS observations of V339 Del 


Date 

2013 

UT* 

(hm) 

MJD 


Exposure 

(s) 

Range 

(nm) 

Resolution 

Sep 23 

6 16 

56,558.261 

(+38 d)^ 

720 

411-686 

90,000 


8 12 

56,558.342 


900 

667-936 

90,000 


10 07 

56,558.423 


3,000 

303-463 

90,000 

Oct 02 

5 02 

56,567.210 

(+47 d)^ 

3,000 

303-463 

60,000 


6 29 

56,567.271 


600 

411-686 

90,000 


7 18 

56,567.305 


900 

667-936 

90,000 

Oct 03 

9 21 

56,568.390 

(+48 d)^ 

3,000 

303-463 

60,000 

Oct 07 

5 05 

56,572.212 

(+52 d)^ 

4,800 

303-463 

60,000 


7 47 

56,572.324 


960 

411-686 

90,000 


8 17 

56,572.346 


1,500 

667-936 

90,000 


* UT is the universal time at the start of an exposure. 
tDays after the optical (V ) maximum (MJD = 56,520.25). 
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